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Abstract 

To detect a warm-hot intergalactic medium associated with the large-scale structure of the universe, 
we observed a quasar behind the Virgo cluster with XMM-Newton. With a net exposure time of 54 ks, we 
marginally detected an OVIII Ka absorption line at 650.9i° g eV in the RGS spectra, with a statistical 
confidence of 96.4%. The observed line center energy is consistent with the redshift of M87, and hence the 
absorber is associated with the Virgo cluster. From the curve of growth, the OVIII column density was 
estimated to be ^ 7 x 10^^ cm~^. In the EPIC spectra, excess emission was found after evaluating the hot 
intracluster medium in the Virgo cluster and various background components. We inspected the ROSAT 
All-Sky Survey map of the diffuse soft X-ray background, and confirmed that the level of the north and 
west regions just outside of the Virgo cluster is consistent with the background model that we used, while 
that of the east side is significantly higher and the enhancement is comparable with the excess emission 
found in the EPIC data. We consider a significant portion of the excess emission to be associated with 
the Virgo cluster, although a possible contribution from the North Polar Spur cannot be excluded. Using 
the column density and the emission measure and assuming an oxygen abundance of 0.1 and an ionization 
fraction of 0.4, we estimate that the mean electron density and the line-of-sight distance of the warm-hot 
gas are ^ 6 x 10^^ cm^'^ and ^ 9 Mpc, respectively. These numbers strongly suggest that we have detected 
a warm-hot intergalactic medium in a filament associated with the Virgo cluster. 

Key words: cosmology: large-scale structure of universe — cosmology: observations — galaxies: 
clusters: individual (Virgo) — galaxies: intergalactic medium — galaxies: quasars: absorption lines. 



1. Introduction 

From several independent measurements — cosmic 
microwave background anisotropy analysis (Spergel et 
al. 2003), primordial deuterium-to- hydrogen ratio (Buries 
et al. 2001), and Lya forest absorption at high redshift 
{z ^ 3) (Ranch et al. 1997) — the cosmic baryon den- 
sity is now converging to J^b ~ Q.Q^hj^, where /lyo = 
Hq/70 kms~^Mpc~^ and Hq is the Hubble constant. 
The local baryon density appears, however, to be far 
lower than these indications (e.g., Fukugita et al. 1998; 
Bristow, Phillipps 1994; Persic, Salucci 1992). Recent 
large-scale galaxy formation simulations predict that most 
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of the baryons (^ 50%) at the present time reside in a 
warm-hot intergalactic medium (WHIM) with tempera- 
tures of 10^-10'' K, and form filamentary structures as- 
sociated with clusters and groups of galaxies (e.g., Cen, 
Ostriker 1999). This warm- hot gas is heated primarily by 
shock heating of gas accreting onto the large-scale struc- 
tm-e (Dave et al. 2001), and its fraction is increasing with 
time. Therefore, detecting those warm-hot baryons is very 
important, not only for settling the 'missing baryon' prob- 
lem, but also for understanding the formation of large- 
scale cosmological structure. 

Considering the elemental abundance and the ioniza- 
tion fraction, O VII and O VIII absorption lines are the 
best tracers for probing the warm-hot gas at 10^-10^ K 
(Perna, Loeb 1998; Fang, Canizares 2000). In fact, these 
lines have been detected with grating spectrometers on- 
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Table 1. Observation mode and net exposure time. 



Instrument 


Mode 


Filter 


Net exposure 


RGS 


Spectrum+Q 




53.8 ks 


EPIC pn 


Full frame 


Thinl 


44.9 ks 


EPIC MOS 


Full frame 


Medium 


56.9 ks 



board Chandra and XMM-Newton, toward PKS 2155- 
304, 3C 273, Mrk 421, and 3C 120 (Nicastro et al. 2002; 
Fang et al. 2002; Rasmussen et al. 2003; Fang et al. 2003; 
McKernan et al. 2003). In most cases, the center energy is 
consistent with a zero redshift, which suggests detection 
of the WHIM in the local group. However, if its redshift is 
zero, distinguishing it from the warm-hot interstellar mat- 
ter in our Galaxy is not trivial (Futamoto et al. 2004) . On 
the other hand, soft X-ray emission in clusters of galaxies 
has been studied (e.g., Kaastra et al. 2003; Finoguenov et 
al. 2003). These authors claimed the detection of a warm- 
hot component with kT ~ 0.2 keV. However, since the 
energy band below 0.5 keV is significantly contaminated 
by emission from the Milky Way halo and the Local Hot 
Bubble, a precise determination of the quantities of the 
warm-hot gas is very difficult with the limited energy res- 
olution of proportional counters and X-ray CCDs. 

In searching for a breakthrough, we decided to probe the 
warm-hot gas associated with nearby clusters of galaxies 
through absorption lines by observing a quasar behind the 
Virgo cluster with the Reflection Grating Spectrometer 
(RGS; den Herder et al. 2001) onboard XMM-Newton 
(Jansen et al. 2001). In this paper, we report on our 
observation and results. We assume /lyg = 1 throughout 
this paper. Errors are quoted at the 90% confidence level, 
unless otherwise stated. 

2. Sample Selection and Observation 

To find an appropriate target, we searched the ROSAT 
All-Sky Survey (RASS) bright source catalogue (Voges et 
al. 1999) and the source catalogue from pointed PSPG 
observations^ for bright quasars behind the Virgo and 
Coma clusters. The target that we selected was LBQS 
1228-^1116. The location of this quasar is (12'^30™54?1, 
11°00'11"), which is 83f3 (or 0.45 Mpc in linear scale) 
south of M87 (figure 1). Because the redshift is z = 0.237 
(Hewett et al. 1995), it is located well behind the Virgo 
cluster. The ROSAT PSPG count rate was 0.12-0.15 cps. 

The observation was performed with XMM-Newton, 
from 2003 July 13 to July 14. The observation ID is 
0145800101. The instrument mode, filter, and the net 
exposure time are summarized in table 1. Note that, al- 
though 100 ks was scheduled, the net observation time 
was only ~ 50 ks due to background flares. 

<http:/ /www.xray.mpe.mpg.de/cgi-bin/rosat/rosat-survey/> 
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Fig. 1. ROSAT PSPG image around M87. The target, 
LBQS 1228+1116, is located 83.3' south of M87. 

3. Data Analysis and Results 
3.1. Absorption Line 

The observation data file (ODF) was divided into two 
files. We performed data reduction using the XMM- 
Newton Science Analysis System (SAS), version 5.4.1, 
with standard parameters. For the RGS, we checked the 
data from a source-free region of CCD 9 for background 
fiares, and accumulated the signal photons only when the 
count rates in this region were less than 0.15 cps for both 
RGSl and 2. The net exposure time was 8.8 ks for the 
first file and 45.0 ks for the second. The background spec- 
tra were produced using the same data sets. All of the 
spectra were binned in four-channel bins (0.042 A), and 
the first and second halves were summed. 

There was no sign of an O VII absorption line at around 
21.6 A, where only the RGSl operates. The 99.7% up- 
per limit of the equivalent width was 2.8 eV at 21.69 A 
(571.6 eV). Therefore, we concentrated on the OVIII re- 
gion. Figure 2 shows the spectra obtained with the 
RGSl and 2, between 18 A and 20 A. There is an 
absorption-line feature at around 19.1 A in both RGSl 
and 2. We fitted the spectra with a power law and a 
negative Gaussian, multiplied by the Galactic absorption 
(iVH = 2.15 X 10^0 cm-2), using XSPEG version 11.2. Since 
the number of photons, especially near the absorption line, 
is small, we used the C-statistic (maximum-likelihood) 
method. When this method is used, a background spec- 
trum cannot be subtracted from the data. Therefore we 
estimated the background level by fitting the background 
spectra of RGSl and 2 simultaneously with a constant 
model, and included it in the fitting model. The results 
are summarized in table 2; the best-fit model is shown 
in figure 2. The line shape was consistent with a narrow 
line, with a 90% upper limit on a width of cr < 5.1 eV. 
Hence, the line width was fixed at 0.1 eV, which is small 
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Fig. 2. Spectra around the Ovill Ko line obtained with 
RGSl and 2. The cross and circle symbols represent the data 
points of RGSl and 2, respectively. The thick and thin solid 
histograms show the best-fit models for RGSl and 2. The 
backgrounds were not subtracted, but were instead simulta- 
neously fitted with a constant model. The background levels 
are shown with dotted curves. For illustrative purposes, chan- 
nels with flickering pixels were masked. 

Table 2. Best-fit parameters of the absorption line. 
Absorption line 



Energy 
Width (cr)* 
Normalization 



650.9 



TTTF" 

'1.9 



eV 



< 5.1 eV (90%) 
(3.9^|g) X 10-6 



photons cm 



Power law 



Index 3.2 
Normalizationt (3.57 ± 0.35) x lO"'' 



C-statistic 
d.o.f 



177.02 
167 



' lor all fits except ior the determination ot the statistical error 

for this parameter, it was fixed at 0.1 eV. 

t In units of photons keV"-*^ cm~^ s~^ at 1 keV. 



compared with the detector resolution (AA 0.06-0.07 A 
or /S.E ~ 2 eV FWHM). When we calculated the errors 
of the absorption-line parameters, we also fixed the pho- 
ton index of the power-law component. The observed line 
center energy and equivalent width were 650.9j:?| eV and 
2.8 respectively. If we assume that this is an O VIII 
absorption line, the energy shift from the rest frame is 
2.7i":9 cV, which corresponds to cz = 1253^369 kms^i. 
This is consistent with that of M87 {cz = 1307 kms~^). 
Note that the systematic error in the absolute energy scale 
of the RGS is 8 mA (rms) or 0.3 eV at 650 eV (den Herder 
et al. 2001), which is much smaller than the statistical er- 
rors. 

We evaluated the statistical significance of the absorp- 
tion line in the following way. The best-fit model without 
an absorption line gives a C-statistic of 181.76. When 
we include an absorption line, but fix the center and the 
width to the energy of the O VIII Ka resonant line with 
cz = 1307 kms^^ and a = 0.1 eV, respectively, we obtain 



the best-fit C-statistic of 177.43. Thus, the improvement, 
AC, is 4.33. The value of AC approximately follows the 
statistic for 1 degree of freedom if the absorption line 
is just a statistical fluctuation (Cash 1979). In order to 
evaluate the statistical confidence of the line detection, 
we produced 10000 simulation spectra using the contin- 
uum model shown in table 2 and the background without 
the absorption line. We then fitted the spectra with and 
without a Gaussian absorption or emission line model. 
We obtained AC > 4.33 for 3.6% of the simulation spec- 
tra. This is consistent with the value expected from the 
distribution (3.7%). We thus conclude that the chance 
probability for detecting such an absorption line at the 
very wavelength corresponding to O VIII Ka at the cluster 
redshift is 3.6%; in other words, the statistical confidence 
of the absorption line is 96.4%. 

The width of the wavelength range where both the 
RGSl and 2 operate is ~ 15.6 A. Because the wavelength 
resolution is ~ 0.06 A, there are about 260 independent 
wavelength bins. The probability for detecting an emis- 
sion/absorption structure of AC > 4.33 at some wave- 
length is almost unity with an expected occurrence of 9.4. 
We searched for such structures in the observed spectra, 
and actually found seven such structures other than that 
at 19.05 A. However, none of them can be identified with 
any atomic emission/absorption line with a reasonable os- 
cillator strength at z = 0, the cluster redshift, or the quasar 
redshift. Thus, these are all considered to be statistical 
fluctuations, and only the 19.05 A feature can be real at 
a statistical confidence of 96.4%. 

We then investigated the systematic errors. First, we 
checked the RGS spectra of bright X-ray sources in the 
archival data to ensure that there is no instrumental fea- 
ture at that wavelength. Then, since the estimation of 
an absorption line could be influenced by the determina- 
tion of the continuum, we examined the dependence of 
the equivalent width on the wavelength region used in the 
fits and the choice of the continuum model. For this pur- 
pose, we changed the width of the wavelength band for 
fits from 1.0 to 2.8 A and added a fifth-order polynomial 
to the continuum. We found that, as long as the width of 
the fitting region is wider than 1.6 A, the variation of the 
equivalent width is smaller than 0.2 eV. When the fitting 
region is too narrow, the continuum level is not well con- 
strained, and is strongly coupled to the equivalent width 
of the line. 

Assuming that the absorption line is not saturated, the 
equivalent width is related to the ion column density in 
the following way: 



EW 



•nhe 



/o; 



(1) 



TOoC 'l-|-z' 

where TVion is the column density of the ion and fos is the 
oscillator strength of the transition (Sarazin 1989). We 
used /os = 0.70 for the O VII resonance absorption line and 
/os = 0.42 for that of O VIII (Verner et al. 1996). Then, 
the O VIII column density was (6.2^4^) x 10^^ cm~^, while 
the 99.7% upper limit on the OVII column density was 
3.7 X cm"2. 
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3.2. Excess Emission 

We also analyzed the EPIC data, and searched for dif- 
fuse emission. Since we are mostly interested in the low- 
energy band, the pn detector is suitable for our analysis. 
Hence, we describe the results obtained from the pn data 
here. We confirmed that the MOS data gave similar re- 
sults, but with larger errors. 

The number of events above 10 keV is regarded as a 
good indicator of the internal background. Thus, we dis- 
carded the data where the rate of pattern events in 
this energy range was larger than 0.7 cps. Then, we ex- 
cluded LBQS 1218-1-1116 and other point sources using 
SAS cdetcct_chain, and accumulated photons from the 
diffuse emission in the 0.3-3 keV band over the field of 
view. To remove any contribution from the internal back- 
ground, we subtracted data taken with the filter wheel in 
the closed position from the spectrum. We then modeled 
it with four components: hot intracluster medium (ICM) 
of the Virgo cluster, cosmic X-ray background (CXB), 
Milky Way halo (MWH), and Local Hot Bubble (LHB). 
For the CXB, MWH, and LHB, we used parameters ob- 
tained by Lumb et al. (2002) as a template. When we 
fixed the temperature of the MWH (0.20 keV), the tem- 
perature of the LHB (0.07 keV), and aU of the CXB pa- 
rameters, the temperature of the ICM determined by the 
fit became 2.0lg'^ keV, which is consistent with that ob- 
tained with ASCA (2.14 ±0.12 keV, Shibata et al. 2001). 
The normalization of the MWH was, however, a fac- 
tor of 2.3-times larger than that obtained by Lumb et 
al. (2002). From an analysis of eight data sets obtained 
with XMM-Ncwton, Lumb et al. (2002) indicated that the 
mean deviation of 0.2 -1 keV flux is about 35% from field 
to field. Therefore, the level of the warm-hot emission 
around LBQS 1228-1-1116 is significantly higher than that 
of the typical Galactic background, even if we consider 
its fluctuation. We flxed the parameters of the MWH, 
LHB, and CXB to those obtained by Lumb et al. (2002), 
added another optically-thin thermal plasma component, 
and flttcd the entire spectrum. Note that the normal- 
ization of the additional component (N) strongly couples 
with the oxygen abundance {A), and A x N is a good 
parameter. Here, we assumed A ~ 0.1. The results are 
summarized in table 3, and the best-flt models are shown 
in flgure 3. Note that the errors given in table 3 do not 
contain the uncertainties of the background components. 
The normalization of the excess emission depends on the 
normalization, abundance, and temperatures of the MWH 
and LHB, while the temperature of the excess emission is 
sensitive to those of the MWH and LHB. 

It is known that the Virgo cluster is located close to the 
North Polar Spur (Loop I) (e.g., Irwin, Sarazin 1996). To 
evaluate its effect, we inspected the RASS 3/4 keV map of 
the diffuse X-ray background (Snowden et al. 1997). The 
upper panel of flgure 4 shows a 25° x 25° image centered 
at LBQS 1228-flll6. The lower panel is a projection of 
the 0.°5 X 50° region shown in the upper panel. The levels 
of the background components (MWH, LHB, CXB) esti- 
mated from the data obtained by Lumb et al. (2002) and 



Table 3. Best-fit parameters of tlie emission spectrum*. 



Intracluster medium (ICM) 
Temperature kT 2T7^^^^TeV 
Normalization iVt (3.0 ± 0.2) x 10"^ 

Abundance A 0.221°;°^ 

Excess emission 
Temperature kT 0.21 ±0.01 kcV 
Normalization (1.5 ± 0.2) x 10"^ 

Abundance A 0.1 (fixed) 

205^36 

d o.f. 146 

" i'he parameters ot the background components were fixed at 
the values obtained by Lumb et al. (2002). Errors do not contain 
systematic errors of the background components, 
t In units of 10~^''/47r(DA(l-l-2))^/ nldV, where Da is the angu- 
lar size distance to the source (cm) and rio is the electron density 
(cm^-^). 
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Fig. 3. Soft X-ray spectrum of the diffuse emission obtained 
with EPIC pn. The aluminum line region was masked. 
Internal background was subtracted using data taken with 
a filter wheel in the closed position. The parameters of the 
cosmic X-ray background (CXB), Milky Way halo (MWH), 
and Local Hot Bubble (LHB) were fixed to those obtained 
by Lumb et al. (2002). The data require another component 
that can be modeled with an optically thin thermal plasma 
with kT = 0.21 keV. 

of the warm-hot emission plus LHB and CXB estimated 
from the present EPIC pn data arc also shown in the panel 
with a dashed line and a solid line, respectively. The level 
of the west side of the Virgo cluster (< —3°) is consistent 
with that of Lumb et al. (2002), while that of the east 
side (~ +4°) is significantly higher and comparable with 
the level of the warm-hot emission plus LHB and CXB at 
LBQS 1228+1116. 

4. Discussion 

Our results suggest the existence of a red-shifted O VIII 
resonance absorption line with a statistical confidence of 
96.4% due to limited statistics. The velocity shift is 



No. ] 



right astenlioii (iimO) 
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■c 
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projetled dislante from LBQS 1 22^+1 1 16 (deg) 



Fig. 4. (Upper panel) 25° X 25° image centered at 
LBQS 1228+1116, obtained from the RASS 3/4 keV map of 
the diffuse X-ray background (Snowden ct al. 1997). The con- 
tours represent 80, 120, and 180 X 10~^ counts arcmin"^. 
(Lower panel) Projection of the 0.°5 X 50° region shown in 
the upper panel. The solid line represents the level of the 
warm-hot emission plus LHB and CXB estimated from the 
present EPIC data, while the dashed line represents the 
background level (MWH, LHB, CXB) estimated from Lumb 
et al. (2002). 

12531 

369 which is consistent with that of M87 

{cz = 1307 kms-i). Thus, we claim detection of the 
WHIM associated with the Virgo cluster. The equivalent 
width of the OVIII Ka line is 2.8^2:0 ^V. This is much 
larger than those reported for WHIM in the local group 
(EW = 0.1-0.4 eV; e.g., Rasmussen et al. 2003; Fang et 
al. 2002). The intrinsic width of the line was not resolved, 
with an upper limit of cr < 5.1 eV (90%). This corresponds 
to a Doppler b parameter of < 3300 kms~^. Because the 
Virgo cluster is thought to be elongated along the line of 
sight by 12-30 Mpc (Yasuda et al. 1997), a velocity differ- 
ence of 800-2100 kms~^ is expected due to cosmological 
expansion. On the other hand, a turbulent velocity of 
a few hundred kms~^ is expected for WHIM from sim- 
ulations (e.g., Ccn. Ostriker 1999). From the curve of 
growth of the O VIII Ka absorption line (see Futamoto et 
al. 2004), the O VIII column density -/Vq vni is estimated to 
be 6.8 X 10^^, 9.6 x 10^^, and 2.4 x 10^^ cm'^, respectively, 
for b = 2100, 800, and 400 kms'^ with EW = 2.8 eV. 
-^oviii ~ 1 X 10^^ cm^^ is close to the maximum column 
density predicted by Fang and Canizares (2000), who in- 



5 

vestigated the column density of the WHIM gas by Monte- 
Carlo simulations. 

On the other hand, no O VII absorption line was de- 
tected. From the upper limit, the ionization fraction ratio 
of OVIII to OVII is constrained to be > 1.7, assuming 
that both of the O VIII Ka and O VII Ka lines are not sat- 
urated. If we further assume coUisional ionization equilib- 
rium, this impHes kT > 0.20 keV (Mazzotta et al. 1998). 
However, photoionization by the cosmic X-ray and UV 
background radiation would increase the ionization frac- 
tion of O VIII if the density is as low as 10~^ cm~'^ (Chen 
et al. 2003). Since we can constrain only the upper limit of 
the density (see next paragraph) , this could be the case for 
the WHIM around the Virgo cluster, and its temperature 
could be lower. 

Excess emission was found in the EPIC spectra over the 
multi-component model representing the emission from 
the Virgo ICM, CXB, LHB, and MWH. The intensity of 
the multi-component model is consistent with the RASS 
diffuse soft X-ray background map in the north and the 
west regions just outside of the Virgo cluster. This 
strongly suggests that the excess emission is associated 
with the Virgo cluster. However, a possible contribution 
by emission from the North Polar Spur cannot be excluded 
completely. We thus take the emission intensity as the 
upper limit for the emission from the WHIM around the 
Virgo cluster. The O VIII column density is related to the 
average electron density, nc, and the line-of-sight length, 
L, of the WHIM as A^oviii = /-AZqUcL, where / is the 
ionization fraction, A is the oxygen abundance relative 
to the solar value, and Zq is the solar abundance of the 
oxygen. On the other hand, if we assume a uniform dis- 
tribution of the WHIM, the emission measure, EM, of the 
plasma is EM ^ n^LS, where S is the geometrical area 
that we observed; and is 7.3 x 10'^ kpc^ for the present 
data. Then, the electron density and the line-of-sight dis- 
tance are constrained to 




V6.2 x 1016 cm-2y ' 

where N is the normalization of the excess emission (see 
table 3), and 0.4 is the ionization fraction for kT ~ 
0.21 kcV. The derived density corresponds to a baryon 
overdensity of i5 ^ 250. The depth of ^ 9 Mpc is much 
larger than the linear dimensions of the Virgo cluster 
in the sky, but is consistent with the elongation of 12- 
30 Mpc suggested by Yasuda, Fukugita, and Okamura 
(1997), which is understood as being a filament of the 
large-scale hierarchical structure along the line of sight. 
Note that the line broadening due to the Hubble flow 
of the 9 Mpc line-of-sight distance corresponds to about 
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1.4 eV. This is much smaller than the 90% upper limit of 
the line width obtained from the present data. However, 
it could be resolved if the statistics were good enough, 
since it is comparable to the energy resolution of the RGS 
(-2 eV FWHM). 

In summary, the present results strongly suggest the 
presence of the WHIM in a filament associated with the 
Virgo cluster at the 96.4% confidence level. To further 
constrain the physical parameters of the WHIM, a precise 
determination of the absorption line profile is an impor- 
tant next step. Future high-resolution spectroscopic ob- 
servations of the emission are also strongly desirable for 
resolving the Galactic and extragalactic components. 

We acknowledge Dr. Yuichirou Ezoe, who helped in 
the analysis of the RASS data. This work was sup- 
ported in part by Grants-in-Aid for Scientific Research 
by MEXT/JSPS (14204017, 12440067, and 15340088). 
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